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abstract  rCanihiua  an  raratma  alia  U nacaaaarr  aaC  lOaniltr  Or  Mack  namkat) 

A new  type  of  instability,  similar  to  thermal  convective  instabilities,  has  been  found  to  occur 
in  self-gravitating  spherical  clouds  uniformily  expanding  or  collapsing  in  a vacuum.  Unlike  the 
Jeans  instability,  it  has  a growth  rate  which  increases  with  zonal  mode  number  I,  so  that  objects 
of  arbitrarily  small  size  are  unstable.  The  modes  are  compressional  and  have  non-vanishing  vorticity. 
The  eigenfunctions  peak  at  the  boundary,  so  the  perturbations  develop  predominantly  into  surface 
inegularities.  The  perturbation  amplitudes  grow  algebraically  at  first,  then  at  late  times  become  — p 
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proportional  to  the  cloud  radius  For  a collapsing  system  in  which  stars  are  produced  by  the 
condensation  of  diffuse  matter  the  theory  predicts  fragmentation  and  enhanced  condensation  in 
the  peripheral  regions.  _ 
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LINEAR  STABILITY  OF  SELF-SIMILAR  FLOW: 

5.  CONVECTIVE  INSTABILITY  IN  BOUNDED 
UNIFORM  SELF-GRAVITATING  SPHERES 

In  this  letter  we  discuss  a type  of  fluid  instability,  the  'convective  gravitational  instability,” 
which  can  occur  in  a nonstationary  self-gravitating  system. 

Unlike  the  Jeans  instability  (see,  e.g.,  Spitzer,  1978),  it  does  not  involve  a critical  mass  or 
size,  and  short,  rather  than  long  wavelengths,  are  the  most  unstable  (Book  and  Bernstein, 
1978).  It  is  also  distinct  from  the  Rayleigh-Taylor  instability,  and  in  a spherical  system  can 
arise  even  when  the  pressure  and  density  gradients  satisfy  pXr)  pXr)  > 0 everywhere.  Rather, 
it  is  closely  related  to  the  convective  instability  found  in  a thermally  stratified  fluid  with  a tem- 
perature profile  such  that  g ■ V T > 0,  where  g is  the  gravitational  acceleration  (see  Landau  and 
Lifshitz,  I9S9).  This  convective  instability  persists  even  in  the  limit  of  vanishing  gravitational 
interaction,  a physical  situation  discussed  by  Book  (1979).  There  it  was  established  that  the  cri- 
terion for  instability  in  an  ideal  polytrope  is  sXr)  <0  [$(r)  is  the  entropy  density].  It  can  be 
shown  (Book  and  Bernstein,  1980)  that  this  same  inequality  governs  the  onset  of  instability 
when  self-gravitation  is  included. 

Here  we  consider  the  convective  gravitational  instability  of  a particular  model  (Staniuko- 
vich,  1949)  which  turns  out  to  be  completely  soluble:  a sphere  of  polytropic  fluid  with  y - 4/3, 
having  uniform  density,  and  undergoing  uniform  self-similar  expansion  (or  contraction)  under 
the  combined  influence  of  inertial,  pressure  and  Newtonian  gravitational  forces.  As  in  related 
stability  problems  of  self-similar  flow  treated  previously  by  Bernstein  and  Book  (1978)  and 
Book  (1978,  1979),  we  employ  Lagrangian  variables.  This  greatly  simplifies  the  analysis,  since 

both  unperturbed  and  perturbed  equation  become  separable. 

Mifiutcript  Mbmitted  March  26,  1979. 
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As  will  be  seen,  this  work  diflers  from  ihe  many  previous  treatments  of  gravitational  in- 


subility  primarily  in  reuining  both  finite  pressure  and  finite  radius  in  the  model.  The  instabil- 
ity, which  is  localized  at  the  periphery  of  the  sphere,  is  driven  by  the  buoyancy  which  fluid  ele- 
ments experience  in  the  combined  inertial  and  gravitational  force  fields.  We  can  therefore 
expect  the  instability  to  give  rise  to  enhanced  mixing  in  the  outer  layers  of  the  system.  This 
result  has  clear  implications  for  star  formation  in  contracting  proto-galaxies  and  possibly  for 
stellar  evolution  as  well. 

We  start  with  the  equations  describing  ideal  hydrodynamic  motion  of  a self-gravitating 


p -t-  V • V p • 0;  (1) 

p V -t-  V p -t-  pV^  — 0;  (2) 

(pp->y  - 0;  (3) 

- 4irCp.  (4) 

Here  p is  mass  density,  v is  velocity,  p is  pressure,  and  0 is  the  gravitational  potential.  Dots 

are  used  to  denote  total  (convective)  derivatives;  y is  the  ratio  of  specific  heats  and  G is  the 

gravitational  consunt.  For  a self-similar  flow  of  the  type  known  as  uniform  or  homogeneous, 

the  radial  position  R at  time  / of  an  arbitrary  fluid  element  whose  initial  position  was  r is  taken 

to  satisfy 

R - r/(/).  (5) 

where  /(O)  1.  Equation  (1)  then  yields 

p(r,f)  - po(r)  (6) 

Then  from  Eq.  (3)  , 

pir.t)  - s(r)  p^  - s(r)  p^  (7) 

where  the  entropy  function  j(r)  is  arbitrary.  In  this  paper  we  will  assume  that  initially  the 

mass  is  distributed  uniformly  inside  a sphere  of  finite  radius  with  density  po.  This  spherical 

cloud  is  surrounded  by  vacuum,  so  that  p(r)  0 — p(r)  for  r > ro  Equation  (4)  can  be 


integrated  to  give 


4(r,t)  - InGp^Vif. 
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Substitution  of  Eqs.  (S)-(8)  in  Eq.  (2)  yields 


f Por^  dr  - 0. 

Equttion  (9)  cm  be  solved  by  separation  of  variables,  yielding 


p(r,/)-f»o(I  -/-Vro^l/-^ 


where  Fo  - s(0)p;. 


The  time  dependence  is  then  determined  by 

(II) 

Po'o  3 

In  analysing  (11),  we  And  that  the  charmter  of  the  solutions  depends  on  he  magnitude  of  y, 
the  relative  magnitude  of  self-gravitational  and  pressure  forces,  and  how  large  / is  in  the  initial 
(/*  I)  stale.  To  simplify  the  analysis,  «ve  set  y ~ 4/3,  obtaining  the  solution  first  found  by 
Staniukovich  (1949).  One  can  dispense  with  either  this  assumption  or  that  of  uniform  density; 
the  analysis  of  the  general  cases,  which  is  somewhat  involved,  will  be  given  elsewhere  (Book 
and  Bernstein,  1980).  It  is  convenient  to  deAne 


lirGpSr^ ' 

and  to  introduce  the  characteristic  time,  deAned  by 


r ^ - 4irGpo/3. 


Equation  (II)  then  becomes 


tV^/-/3-1.  (14) 

The  Arst  integral  of  (14)  can  now  be  written 

y (/*-/o)  - (/5  - 1)  (1 -r').  (IS) 

where  /o  is  the  value  of  /(/)  at  t - 0.  Following  Sedov  (1959),  we  distinguish  three  types  of 
motion.  In  type  I motions,  / 0 for  all  /,  0 < / < oo.  Trsjectories  of  Type  I for  which 

/ *-  0 as  /— *00  become  possible  if  B < 1: 


I 
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-/3)r'.  (16)  I 

I 1 

When  there  is  a turning  point,  we  take  the  zero  of  time  there: 

1 

i 

(17)  ^ 

* I 

which  yields  solutions  of  Type  II  for  /3  < 1 and  Type  III  Tor  /3  > I.  ^ 


i 


[ 


We  return  to  Eqs.  (l)-(4)  and  make  use  of  the  results  obtained  by  Bernstein  and  Book 
(1978)  and  Book  (1978,  1979)  in  deriving  linearized  equations  for  the  evolution  of  a small  per- 
turbation about  the  basic  state  solutions.  The  perturbed  displacement  ( satisfies  the  linearized 
form  of  Eq.  (2), 

p'£  + Pi  R ™ — V nPi  + (V n^)  ' V tP  — pV n^i  + p(V  — Pi  V n^,  (18) 

where  first-order  quantities  are  distinguished  by  the  subscript  1.  Here  the  gradient  operator  Vk 

is  defined  with  respect  to  the  unperturbed  coordinate  R.  The  perturbed  density  is 

Pi--p^Ri  (19) 

and  the  perturbed  pressure  is 

Pi--ypV|,|.  (20) 

Using  Eqs.  (19)  and  (20)  and  the  unperturbed  forms  of  Eqs.  (2)  and  (4),  we  rewrite  Eq.  (18) 

as 


! 


i 


I 


j 


1 


p£+  (Vnp)  V,  yV,  (p  V,  •^)  -t-  p(V,£)  • R + pV,<6i  - 0.  (21) 

We  find  an  equation  for  by  substituting  Eq.  (8)  in  Eq.  (4)  and  collecting  first  order 
terms: 

Vjl^i  — Vj  • “ 4irCp|.  (22) 

At  this  point  it  is  convenient  to  use  Eq.  (S)  to  rewrite  the  equations  in  terms  of  r instead  of  R ' ' 

and  then  to  introduce  the  notation  V = p 

We  can  eliminate  ^|  from  Eq.  (21)  either  by  operating  with  the  curl  or  by  taking  a diver-  ^ 

gence  and  substituting  from  Eq.  (22).  The  results  are  given  by  ! 
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(1  - r^a  + 2(1  + /J)<r  +/3« 


and 


/J  ii  - O - 1)  « + pill  ■ Vtr  + rr:VVa  - V). 

where  r is  expressed  in  units  of  tq  and  t is  units  of  r,  and  we  have  defined 


«r  - V 

«i»-r-VxVxf. 


We  seek  separable  solutions  of  the  form 


(Tit.i)  - S(r)Y,„{9.4>)  TU) 
w(r.i)  - W(r)Y,j0.4»)  TO). 


Eliminating  W,  we  find 


r-  (/3-  1)  fiT. 


and 


/■-^jlrMCl  - r^lSir-  /(/  + 1)  (1  - r^)s]  + ATS  - 0. 


Here 


K - 


2/3(1 


and  M is  a separation  constant  which  determines  the  time  behavior.  Let 

z(r)  - /-(I  - r^)S. 


Then  Eq.  (30)  can  be  written 


z"  + 


K 1(1  + 1) 


z -0. 


(23) 


(24) 


(25) 

(26) 


(27) 

(28) 


(29) 


(30) 


(31) 


(32) 


(33) 


1 -r^  ^ 

Since  Pi  must  be  finite  at  r - 0 and  vanish  at  r - 1,  the  boundary  conditions  appropriate  to 
Eq.  (33)  are  z(0)  - z(l)  - 0. 


If  one  multiplies  Eq.  (33)  by  z,  there  results  after  an  integration  by  parts 


< 4 

1 


S 


BERNSTEIN  AND  B(X)K 

K -1 > 0.  (34) 

f'  drz^ 

Jo  1 

Note  that  if  we  difTerentiate  Eq.  (31)  with  respect  to  m for  Axed  /,  it  follows  that 

Anticipating  a result,  proved  below,  that  the  rate  of  growth  of  an  unstable  mode  satisfying 
Eq.  (29)  is  greater,  the  nearer  /x  (/3  - 1)  approaches  +oa,  we  see  that  for  a given  / maximum 
growth  is  obtained  for  the  mode  with  the  smallest  value  of  K. 

Regularity  of  z at  /•  — 0 implies  via  Eq.  (33)  that  for  small  r,  z(r)  — r'^'.  Since  the 
equation  is  invariant  under  r - r,  one  can  seek  a solution  in  the  form 


where  the  {c„)  are  constant.  The  recursion  relation  is  readily  seen  to  be 

_ A — Anin  + / + 1/2) 

“ 4(»  + 1)  (n  + / + 3/2)  "■ 

If  K is  chosen  such  that  the  series  (36)  terminates,  a solution  for  which  z(l)  — 0 can  be 
obtained  (since  in  Eq.  (33)  the  second  term  would  otherwise  diverge].  The  condition  for  this 
to  occur  is 

K •AMN  + I + 1/2),  (38) 

N — 1,2,...,  The  smallest  eigenvalue  is  clearly  A — 4/  + 6,  corresponding  to  N * 1.  The 


associated  eigenfunction  is 


or  equivalently,  S(r)  - r'. 


z(r)  - r'*'  (1  -r^. 


Taking  A — 4/  + 6 in  (31)  and  solving  for  fi,  we  get 
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M-  I 


3 4ip 

3 , Aip 

^ . HI  +l)P^ 

2 3(/3-l)  * 

2 3(/3-l) 

* W-l)’ 

(40) 


With  the  appropriate  sign,  (m  - D (/3  - I)  > 0 for  all  / > 0.  This  is  the  condition  for  instabil- 


itity,  as  will  now  be  shown. 


The  usual  definition  of  instability  is  inadequate  when  applied  to  nonsteady  states,  since 
the  time  dependence  Tit)  of  the  perturbations  is  in  general  not  exponential.  It  is  appropriate 
to  call  a mode  stable  (unstable)  if  the  ratio  dit)  - |r(r)//(r)|  of  the  perturbation  amplitude 
to  that  of  the  basic  state  vanishes  (diverges)  as  / — • oo.  Note  that  this  definition  is  useful  both 
when  / “•  0 and  when  / — • oo  for  large  i.  Moreover,  it  does  not  preclude  the  possibility  that, 
e.g.,  in  an  ’unstable*  system  Qit)  decreases  for  a finite  time.  Because  the  period  during  which 
/ is  effectively  nonzero  is  finite,  Q (<»)  may  be  finite.  If  this  happens,  the  pertinent  question 
becomes  that  of  finding  how  much  relative  amplification  occurs,  i.e.,  of  evaluating  d (<»)/fi(0). 

If  ip  — 1)  t*  » 1 in  Eq.  (29),  one  solution  grows  and  the  other  damps  approximately 
exponentially  for  / = I,  experiencing  roughly  e-foldings.  Since  / = I holds  only  for 
times  |r|  < 1/2,  this  is  a very  crude  measure  of  the  real  time  amplification.  Fortunately,  it  is 
possible  to  solve  Eq.  (29)  exactly.  This  is  easily  accomplished  by  using  / instead  of  t as  the 
independent  variable. 

For  Type  I motions  Eq.  (29)  becomes  by  virtue  of  Eq.  (16) 

2/'-0-/-^+Mr-O.  (41) 

o/'  a/ 

Hence  Tit)  is  a linear  combination  of  the  two  solutions 

TJi)  - (42) 

where  A - (9  - 8ft) For  /»  < 9/8,  A is  real  and  grows  more  rapidly  than  r_  with 

increasing  /.  If  ft  < 1 , 


(43) 
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diverges  is  / — «».  while 

vanishes.  For  outward  motion,  evidently  perturbations  with  a component  proportional  to  7*^ 
blow  up  at  late  times  faster  than  the  unperturbed  radius,  and  consequently  satisfy  our  condition 
for  instability.  If  m > 1.  on  the  other  hand,  the  perturbations  are  stable.  Thus  the  compressi- 
ble modes  given  by  Eq.  (29)  with  n determined  from  Eq.  (40)  can  be  stable  or  unstable, 
depending  on  which  branch  is  taken.  For  inward  motions,  (2_  grows  more  rapidly  than 
Both  diverge  when  n > 9/8,  albeit  slowly.  This  divergence  has  been  related  by  Book  (1978)  to 
the  fact  that  y < S/3.  The  explanation  is  essentially  geometrical,  and  the  limiting  behavior  can 
be  derived  by  invoking  conservation  of  the  action  of  a standing  wave  on  the  collapsing  cloud. 
When  M " I.  f When  m < K T+  is  well  behaved  as  / — 0,  but  T.  is  unstable.  A dis- 
turbance initialized  at  any  instant  of  time  with  spatial  dependence  ~ r**'  blows  up  as  the 

cloud  collapses. 


For  Type  11  motions  we  find  that  as  a function  of  jt  — 1 - /■',  T satisfies  the  hyper- 
geometric equation 

xd-xly”-!-  Ic  - (o  -I-  b -h  Dxly'-  aby  - 0 (45) 

with 

"I  - -J-  (3  ± A);  (46a,b) 

c - y-  (46c) 

Writing  the  two  independent  solutions 


?(/)  - jF, 


y(3-l- A).  ■J-(3- A);  y;l  -/  ' 


§(/)  - I2(/-' - 7(5  + A),  7(5-A);|;1  ~ 

4 4 2 


(47) 

(48) 
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SO  that  in  terms  of  the  initial  values. 


•(3+At/4 


When  n > 9/8,  A is  imaginary.  The  second  term  in  each  of  the  above  equations  is  the 


complex  conjugate  of  the  first,  and  all  terms  oscillate  for  small  /,  diverging  as  / When 


H < 9/8,  A is  real  and  the  second  term  in  both  equations  dominates  as  / — > 0.  For  — 1 


?(/)  “ /(r).  Note  that  the  dependence  on  /in  Eqs.  (52)-(53)  is  asymptotically  identical  with 


that  of  the  exact  solutions  given  in  Eq.  (42)  for  Type  I motions. 


For  Type  III  motions  the  solutions  J and  g are  still  valid,  with  a,  b and  c as  defined  in 
Eq.  (46).  We  are  interested  now  in  the  stability  of  perturbations  on  a cloud  expanding  outward 
(/—»<»).  The  appropriate  asymptotic  forms  are 


BERNSTEIN  AND  BOOK 

The  coefficients  /"*“*  and  G'°°*  increase  monotically  with  /x.  Although  the  amplification  for 
both  Eq.  (54)  and  Eq.  (55)  remains  finite,  it  can  be  made  arbitrarily  large  as  ^ • oo.  We  thus 
see  that  the  modes  undergo  a growth  which  can  be  catastrophically  large. 


If  we  examine  the  reverse  trajectory  (the  cloud  collapsing  from  infinity,  turning  at  / » 1, 
and  rebounding)  we  find  a clearcut  instability.  The  solution  of  Eq.  (29)  having  unit  amplitude 


as  / — > - oo  IS 


K(,)-f  5^, (56) 


5 + A 


5 - A 


/”*  1-  Hence  the  relative  amplification  3 (0)/3(—  °°)  is  infinite.  If  the  perturbation  is  ini- 
tialized when  the  cloud  is  at  some  large  but  finile  radius  /(—/(»),  the  relative  amplification  at 
r « 0 of  mode  /is  proportional  to  /(-ro)  //'"*  (/),  which  increases  with  both  /and  /3. 


By  choosing  an  appropriate  linear  combination  of  the  two  independent  solutions  of 
Eq.  (29),  a solution  can  be  found  for  motions  of  Type  I,  II  or  III  which  decreases  initially,  then 
grows.  This  justifies  our  choice  of  the  asymptotic  relative  amplification  as  a criterion  for  instabil- 
ity. 


1 


In  conclusion,  we  note  that  the  preconditions  for  instability,  viz.,  that  the  outermost 
layers  of  the  system  be  cooler  than  the  core  (str)  < 0],  is  not  at  all  implausible.  In  astrophy- 
sical  contexts,  it  might  result  from  surface  radiation  or  from  nucleosynthesis  in  the  dense  core. 

Even  when  p « I,  short-wavelength  surface  perturbations  would  grow  rapidly,  leading  to 
mixing,  condensation  of  small  masses  and  development  of  fine-scale  structure  there. 

j 
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